Continuous 3d Freezing Transition in Layered Superconductors 
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We use Ginzburg-Landau theory to study the H C 2 transition in layered superconductors with field 
parallel to the layers, finding a continuous 3d freezing transition to a triangular vortex super-solid 
in the three-dimensional XY universality class. If screening effects are neglected, off-diagonal-long- 
range-order survives only for d > di c = 5/2. The partial breaking of the lowest Landau level 
degeneracy induced by layering leads to a local selection of a triangular lattice structure, in contrast 
to the global free energy minimization in, e.g. Abrikosov's calculation. Our results are relevant to 
artificially layered superconductors and to strongly anisotropic high T c materials. 

PACS: 64.60.Fr, 74.20. D 



The development and application of the renormaliza- 
tion group in the 1970's and early 1980's greatly im- 
proved the understanding of phase transitions. Despite 
many dramatic successes resulting from this creative ex- 
plosion, the normal to superconducting transition in fi- 
nite field (NSFF) , one of the earliest formulated examples 
of critical behavior, is still not understood. 

Any description of the NSFF transition begins with 
the Ginzburg-Landau (GL) free energy 
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where ro ~ (T — T c ) /T c , (f>o = hc/2e is the flux quantum, 
fio is magnetic permeability of the normal metal and is 
close to unity, and we have for simplicity ignored the tech- 
nically unimportant effective mass anisotropy in the first 
term. Eq^ respects several symmetries which are broken 
within mean-field theory in the mixed phase: translation 
and rotational invariance, as well as the C/(l) gauge sym- 
metry corresponding to off-diagonal long-range-order 
(ODLRO) and the choice of phase of the pair wave- 
function Probably the most consistent definition of 
a 'gauge-invariant phase' was suggested by Moore jj], 
and is completely specified by the positions of vortices, 
or zeros of "J. A hidden third symmetry was introduced 
in Ref. m as a (different!) U(l) phase invariance of the 
dual "boson" (dis-)order parameter tp. The expectation 
value ("0) ^ in the normal phase. 

The mean-field description of the NSFF transition was 
discovered by Abrikosov 0]. For fields smaller than 
H = H C 2(T), both 'i'(r) and B become periodic, with 
an amplitude (of \& and of the periodic modulation of 
B) increasing continuously from zero as H or T is re- 
duced. While a consistent (i.e. non-mean field) theory 
of the phase transition is still lacking (but see below), 
fluctuation effects are well understood within the lattice 
phase. As discussed in Ref. , phonon fluctuations of the 



vortex lattice destroy ODLRO for d < 4. However, the 
mixed phase can persist with broken translational (and 
rotational) invariance. In principle, two types of lattice 
phases can exist: the ordinary solid, with (-0) = 0, and 
the supersolid, with (ip) ^ |||. 

The fluctuation corrected critical behavior has been 
studied in two ways. The more phenomenological ap- 
proach attempts to construct a theory in terms of vor- 
tices, the low energy degrees of freedom of the mixed 
phase. In this picture, the NS transition is accomplished 
by a first order freezing of vortex lines from a liquid 
(the normal phase) into a solid (the vortex lattice) ||[| . 
Somewhat more microscopic treatments begin with the 
Ginzburg-Landau theory and attempt to capture fluctu- 
ation effects near the transition directly in terms of 'L. 
As first pointed out in Ref. 0, the infinite degeneracy 
of the lowest Landau level (LLL) requires a functional 
renormalization group treatment of the critical behavior 
in such models in high dimensions |^||. Unfortunately, 
the full connection between this approach and the more 
physical picture of vortex lattice melting remains unclear. 

In this letter, we consider the case of a layered su- 
perconductor with magnetic field oriented parallel to the 
layers. Because of the partial breaking of LLL degeneracy 
imposed by this intrinsic periodicity, we are able to carry 
the microscopic treatment much further than in previous 
fluctuation studies. This allows a direct connection to 
be made between the Ginzburg-Landau and vortex the- 
ories. It also provides an analytic understanding of the 
selection of a triangular vortex lattice: in all the previous 
studies to date the breaking of degeneracy and ordering 
into a particular lattice has only been captured by an 
explicit numerical comparison of free energies. 

Our most important result is the demonstration of a 
continuous mechanism for three dimensional vortex freez- 
ing in layered superconductors in parallel fields. 

The effect of the Cu02 layers (taken normal to the y 
axis) can be modeled using Eq.[l] and a spatially modu- 
lated local temperature ro = to — acos (2-Ky/d), where 
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d is the layer spacing. For simplicity we first ignore 
screening (or equivalently consider length scales below 
the London penetration length A). With the magnetic 
field B — Bz along z-axis (or more generally transverse 
to the xy plane, in d = d± + 2-dimcnsions), it is conve- 
nient to work in the Landau gauge with A = Byx. The 
quadratic part of the free energy is simplified by expand- 
ing tp in the eigenfunctions of the corresponding equation 
of motion. In this gauge, the eigenmodes are of the form 

* = e ika,x+ik ± -r ±( i with 
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E(k x ,k ± ) = E(k x ) + k\ and £ = y^/2^B. For a = 
the "Hamiltonian" is diagonalized by Landau level (LL) 
eigenfunctions, with the usual LL degeneracy of E with 
respect to k x . The effect of the periodic potential is 
to break the continuous translational symmetry along y 
down to a discrete group, which implies E(k x + k n , k±) = 
E(k x , kx), where k n — dn/£ 2 and n is an arbitrary inte- 
ger. The original degeneracy N$ — L x Ly/2n£ 2 is thereby 
reduced to Nj — L y /d. Each LL is broadened into peri- 
odic bands with width proportional to a. For our further 
approximations to hold, we consider a strong field limit 
which consists of a number of constraints. First, we take 
h2eB/Mc >> a (see below), which implies that LLs re- 
main narrow on the scale of their separation, and allows a 
clear separation of LLL physics. Secondly, we require the 
periodic potential to vary smoothly on the scale of mag- 
netic length, i.e. a << d 2 /£ A . In this limit, the functions 
4>(k x ,kj_,y) are Gaussian, 
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where J k fe = / dk x d d± kx/{2n) (d±+1) . In the strong 
field limit the only effect of the periodic potential is to 
introduce periodic k x dispersion, corresponding to the en- 
ergy cost of moving the center of the Gaussian along the 
cosine potential. The quadratic part of the free energy 
becomes 
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where E(k x , k±) = k\ + to + a cos (2ir£ 2 k x /d) and a — 
ae -(iri/d) 2 _ is convenient to work in the reduced 
zone scheme where k x = k + dn/£ 2 , and define new 
fields 4>n(k,kx) = <f>(k + k n ,k±). The nonintegral part 
k <C d/£ 2 describes small deviations from the local mini- 
mum of the cosine potential. Rewriting Fo in terms of the 
4> n yields a dj_+l-dimensional A^-component field theory, 
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where 7 = (a/2)(2ir£ 2 / 'd) 2 and t = t Q + £~ 2 - a. The ro- 
tational invariance in the iV^-dimensional n-isospin space 
of F Q is of course broken by the the full GL free energy, 
to which we now turn. 

The quartic term is readily re-expressed in terms of the 
</>„(&:, kj_) fields using Eq.||. The result is 

Fiat =g X! / Kl( k l)<l ) n2(k2)(f>n3(k3)(f>n4(k4) 

xV[{n l },{k i }](2n)^^ d - 1 \k 1 -k 2 + ks-k 4 ), (6) 
where 
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Here n\ = n/2 + k/2,n 2 = n/2 + l/2,n 3 = n/2-fc/2, and 
rM = n/2 — 1/2, and we neglected small 0(k) corrections 
in the exponential in Eq.f^. 

Infrared divergences emerging below d uc = 5 can be 
studied using renormalization group (RG) methods in 
d = 5 — e dimensions. We use the momentum shell 
method, in which short wavelength modes in a shell 
A/b < p < A are progressively integrated out to yield 
a renormalized free energy. The RG is best formulated 
in terms of the coupling matrix uu,i — V[ni,n2,Ti3,ni\. 
The matrix Uk,i vanishes when k, I are not simultaneously 
odd or even. The discrete translational invariance along 
y enforces the independence of Uk,l on n. 

We perform an infinitesimal RG transformation with 



and rescale r — > br, 



b^- d )/ 2 (j) n to obtain 



dtUk,i = £u kji — 2mu kjS u S} i — 8u SjS _ k +i u s _ 3k / 2 +i/2,s-k/2-i/2 

(8) 

in the thermodynamic (Nd — > 00) limit. For later pur- 
poses (see below) have considered an extended U(l) x 
0(m) theory with the generalized interaction fi nt oc 
d(fci)^ 2 ( fc 2)C3( fc 3)^ 4 (fc4), with i,j = l,...,m. 

We have succeeded in finding only one exact (but un- 
stable) solution to the Eq.||with u* kl = e/(2m + 8)SkfiSk,i, 
corresponding to simultaneous but independent XY tran- 
sitions within each layer. For m — > 00, the second nonlin- 
earity of Eq.|| becomes negligible and the resulting ma- 
trix equation is easily solved. The large m fixed point 
solutions are then any matrix with an arbitrary combi- 
nation of and e/2m eigenvalues. The the most stable 
solution is proportional to 5^,1, corresponding to an ar- 
bitrarily long-range interaction cx J^ n ,k $n 4h,<t>n+kK.+k- 
However, the form of the initial condition for u° k , = uj^ui, 
Eq.|?], dictated by the LLL physics, insures that the phys- 
ical fixed point solution is of the same short range form, 
but with a single nonzero eigenvalue flowing to e/2m. At 
this fixed point standard methods give the correlation 
length exponent v = 1/(2 — e/m), which diverges (for 
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to = 1) in d = 3, naively suggesting 3 as the lower critical 
dimension (d; c ). Although this result is in accord with di- 
mensional reduction arguments of d — > d±_ = d—1, we will 
show below that a more careful analysis leads to di c = 2.5 
for the unscreened case, invalidating dimensional reduc- 
tion. The breakdown of dimensional reduction may also 
be important for the NSFF transition in an isotropic su- 
perconductor. Unfortunately, the 1/to corrections to the 
to — > oo fixed points appear to be singular, and we have 
not succeeded in expanding toward to = 1. 

Although we have not been able to solve above RG 
equations, if the quartic interaction remains short-ranged 
(as is physically plausible from the existence of the mag- 
netic length scale I), much insight into this freezing tran- 
sition can be gained from the mft and from studying the 
fluctuations in the low T lattice phase. The short-range 
form of V [{rii}, {fci}] provides the key simplification rel- 
ative to the problem without the layered potential. 

Eq.^| contains interactions between all of the {(f> n } at 
all momenta. We expect the minimal free energy to oc- 
cur for k = 0, i.e. uniform </>„; higher order terms in a 
gradient expansion are also less relevant in the RG sense. 
Because increasing n% — n.3 and n2 — n.4 results in expo- 
nential suppression, it is sufficient to keep only the largest 
few terms in the interaction free energy 
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with w — v = 0. Spe- 
and indicate that 



where u = 2 5 /v / 2^F, w = 4g/V2nPe- d ' 2l and v = 
4g/v27rPe~ d I s - . Since w,v <§; u, there will be a broad 
crossover regime described by Eq.| 
cializing to d = d± + 2 = 3, Eqs. 
in this limit the individual <p n are decoupled and com- 
prise two-dimensional XY order parameters in each layer. 
These develop significant amplitude below the mean-field 
transition at r = 0, where the system shows the usual 
broad specific heat maximum associated with the 2d XY 
model. Physically, the large field limit leads to a strong 
variation of the phase of iji, destroying the Josephson 
coupling between neighboring superconducting layers. 

To study the behavior for r < 0, we ignore (massive) 
amplitude fluctuations and let <f> n = -J\r\jue l9n . Up to 
a constant, the free energy becomes 



F 



£ fdxdyiK x (dJ n ) 2 + K z {dJ n f 
+ucos(0 n+ i+0 n _i-20 n )l, (10) 



where K x = jK z = j\t\/u and v = t 2 v/u 2 . For 
5 = 0, the system undergoes simultaneous Kosterlitz- 
Thouless (KT) transitions within each layer (each n) to 



decoupled quasi-long-range-ordered (QLRO) states for 
K = yJK x K y > Kkt ~ 4/7r. Correlation between the 
phases in the layers is established for v > 0. A simple cal- 
culation using standard methods for the 2d sine-Gordon 
model demonstrates that v is indeed a relevant pertur- 
bation for any K > K* = 3/(47r) < Kkt, and there- 
fore preempts the KT transition p2[ . Instead, a true 
three-dimensional transition occurs at a critical stiffness 
K c < -Kkt (K c in fact approaches Kkt as v — > 0). 

For K > K c , the ordered state is determined by 
the minimization of Eq.|l^. Because v > 0, the ar- 
gument of the cosine should be an odd multiple of tt, 
and a convenient choice of symmetry breaking is 9 n ^ = 
0, 0, tt, tt, 0, 0, • ■ • on successive layers. Examination of the 
zeros of l^l 2 , EqJ| shows that this low energy state is a 
triangular vortex lattice with lattice vectors 2ir£ 2 x/d and 
ir£ 2 ic/d + dy [fl3"|| . The simple determination of the lat- 
tice structure from the local v phase interaction strongly 
contrasts with the situation in the conventional unlay- 
ered GL theory, in which the choice of lattice relies upon 
a complex non-local minimization in the full LLL basis. 

Fluctuations within the ordered phase are described 
by the phase shifts 9 n = 6 n — n fl, which are governed 
by Eq.[l0] with the sign of v changed. To study the long- 
wavelength properties of these phase variations, the co- 
sine may be replaced by a mass, leading to 



+ Kz (d z 9) 2 + K v {dl6f\, (11) 



where k x = K x /d, n z = K z /d, n y — 2m 2 d, and the 
thermally renormalized mass to 2 = v ^K/{e,iTK-e,) _ 
straightforward calculation shows that the phase fluctua- 
tions (\0\ 2 ) are finite in three dimensions, with di c = 2.5, 
showing that GL theory (without screening) displays true 
ODLRO, in contrast to the naive dimensional reduction 
conclusion of c?; c = 3. This free energy occurs in the 
study of liquid crystals H] and is known as m=l, XY 
Lifshitz model (m=l, for one soft k component). The 
model can be rewritten as the zero field gauge theory 



d 3 r U ± \V_ 



ay I 



Va 



7. 1 2 



(12) 



where dy is a "fake" gauge field with only a y-component 
and we have rescaled x and z coordinates. Ref. |Q con- 
cludes that the approach to high temperatures can pro- 
ceed in two stages, via an intermediate "orientationally 
ordered" state. This phase lacks ODLRO, but retains 
order in d y 8 ~ u, and hence retains broken translational 
invariance. It is also known (l4|] that the free energy, 
Eq.[l2|, is self-dual, which implies that the intermediate 
state has long-range correlations in the disorder parame- 
ter ip - i.e. it is a super-solid. We see no a priori reason 
why such a topologically ordered phase may not survive 
below die, e.g. in two dimensions. 

Without screening, the solid to super-solid transition 
is therefore in the 3d, to = 1, XY Lifshitz universality 
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class. By analogy with the 3d normal-to-superconducting 
zero-field transition in the strongly type II limit Jl5| , we 
expect a continuous transition. 

What is the universality class of the normal to super- 
solid transition? As described in Appendix B of Ref. [^6) , 
the only broken symmetry associated with the transition 
is translations along the x axis. A Landau formulation in 
terms of the complex order parameter p± = (^)^e ldx l l ) 
implies the critical behavior is of the 3d XY type (c.f. 
also Rcf. 14 1). The phase of p% is proportional to the 



phonon coordinate u. Note that electromagnetic screen- 
ing does not play an important role here, because long- 
range interactions are screened anyway by the fluid of 
vacancies/interstitials already present in the super-solid. 

The above considerations apply when magnetic vari- 
ations may be neglected. Gauge field fluctuations can 
be included (in the ordered phase) simply by the mini- 



mal coupling prescription 8^9 



d„9 



where a = 



27r(A — Byx)/(/)Q and </>o = hc/2e is the superconducting 
flux quantum. The gauge-invariant free energy is thus 



-Fgaugc = J d 3 rl^K x \d x 8 ~ a x \ 2 + n z \d z 9 — a z 



-K y \d y (d y e - a y )\ 2 



32tt 3 
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Note that, despite its unusual form, the K y term is in- 
variant both under gauge transformations and the global 
shift (c.f. Eq.|TT|) 9 — > 9 + Xy. Eqjl^also encodes a trans- 
verse Meissner effect: minimizing Fgauge with respect to 
9 for a = Byzic gives a divergent free energy density for 
an additional field along the y axis. Perturbing fields 
along x and z are not screened. 

We now choose a gauge with 8 = to determine the 
true low-energy structure of the theory. The result is 
a modification of the usual Anderson-Higgs mechanism 
which occurs in the Meissner phase: the a x ,a z compo- 
nents of the gauge field become massive, but a y continues 
to describe a gapless mode. This remarkable mechanism, 
allowed by the unusual k v term in Eq.[L3], has a direct 
physical interpretation. Once gauge fluctuations are in- 
cluded, only the zeros of \& are gauge-invariant objects, 
and it is their displacements u which are the physical 
degrees of freedom. Because of the potential in the y 
direction, only u x = u is gapless. It is straightforward 



to derive the relation u — —£ (d y 9 
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9 = gauge. Integrating out the massive a x and a z fields 
then gives the elastic free energy 



F et = - I d 3 r \ c ± \V ± u\ 2 + c y \d y u\ 



-c 4 \K x 
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where cj_ 
(d x ,d z ). 



^/(32^ 3 £ 4 ) 



2k v /1 4 1 and V_l 



0q/(5127t 6 £ 4 ) may be important experimentally due to 
the extremely small value of c y in the strongly layered 
limit. We observe that with screening the vortex lattice 
survives down to two dimensions. 

Finally, we note that in an isotropic system the low- 
energy degrees of freedom are also a discrete subset k n 
of the LLL modes. Even without layering, rotational in- 
variance in the plane perpendicular to field and the LLL 
physics unambiguously determine k n = n^/2~Tr/l. One 
may speculate that fluctuations within such a discrete 
set of modes could support a continuous freezing mecha- 
nism within the Ginzburg-Landau theory of an isotropic 
superconductor as well. 
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